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Abstract. In the solution of the superintegrable chiral Potts model special polyno- 
mials related to the representation theory of the Onsager algebra play a central role. 
We derive approximate analytic formulae for the zeros of particular polynomials which 
determine sets of low- lying energy eigenvalues of the chiral Potts quantum chain. These 
formulae allow the analytic calculation of the leading finite-size corrections to the en- 
ergy eigenvalues without resorting to a numerical determination of the zeros. 



1 Introduction 

The chiral Z/v-symmetrical Potts model, apart from showing a rich phase struc- 
ture, has attracted considerable interest because for appropriate parameter choi- 
ces it has very special integrability properties. In the "superintegrable" case it 
is integrable both because it provides a representation of Onsager's algebra (in 
this respect being a natural generalization of the Ising model) , and also because 
its Boltzmann weights satisfy a new type of Yang-Baxter equation. Both prop- 
erties guarantee the existence of an infinite set of commuting charges. In order 
to calculate the spectrum, functional relations are needed. Such relations had 
first been conjectured by Albertini et al. (1989), and later derived by exploring 
relations to the six- vertex model at q — e 27Tl / N (Bazhanov and Stroganov 1990, 
Baxter et al. 1990). The formulae which give the spectrum in the superintegrable 
case have been given by Albertini et al. (1989) for Z3 and for general Zn by 
Baxter (1994). Little is known about correlation functions. For the general in- 
tegrable case (not satisfying the Onsager algebra), the free energy and interface 
tensions have been obtained in the thermodynamic limit (O'Rourke and Baxter 
1996). It is still a challenge to find an analytic derivation of the order parameter 
(Baxter 1998). Recently, general integrable boundary conditions for the chiral 
Potts model have been obtained by Zhou (1998). 

In the following we concentrate on the solution of the superintegrable case, 
in which special polynomials related to the representation theory of the Onsager 
algebra play a central role. The energy eigenvalues are determined in terms of 
the zeros of these polynomials. For the lowest states the polynomials are of order 
rriE = [((A — 1)L — Q)/N]. L denotes the number of sites of the chain and Q 
the .Z/v-charge sector, [x] stands for the integer part of a rational number x. 
In the limit L — > 00 the sum over the terms depending on the zeros can be 
expressed as a contour integral without explicitly calculating the zeros. This 
method has been used to obtain the ground state energy and excited levels 
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in the thermodynamic limit. For the gapless regime finite-size corrections have 
been studied by Albertini and McCoy (1990) in order to obtain information 
on correlation functions through arguments of conformal invariance. Interesting 
effects due to the broken parity invariance have emerged. For finite L one has to 
resort to numerical methods. 

In these lecture notes, after reviewing main features of the model and its 
solution, we study the zeros of the above mentioned polynomials. We derive 
quite accurate approximate analytic expressions for their location at finite L. 

2 Integrability of the chiral Potts model 

2.1 The integrable model 

For many discussions of the chiral Potts quantum chain, it turns out to be 
sufficient to consider the following Z^y-symmetric hamiltonians which depend 
on three real parameters A, <p, <p: 

H = -YY J-_ f e M2l~N)/N Z l + A ^4>(2l-N)/N X l X N-1\ 

^ ^ sm(nl/N) \ 3 / 

The first sum is over the sites j. At each site j there are operators Xj and Zj 
acting in a vector space C N and satisfying 

ZiXj = Xj Z % ■ Zf = Af = 1; w = e 2m ' N . 

Often it is useful to represent the operators Xj and Zj by [Xj)^ m — S^ m+ i 
(mod N) and Zj = rfjaj(l,o),ii) 2 , . . . , uj N ~ 1 )j. H commutes with the Zjy-charge 
operator Q — Ylj =1 Zj. We write the eigenvalues of Q as where Q = 
0, 1, . . . , N — 1. A is the inverse temperature. We shall consider the periodic case 
Xn+i = Xi, although, as shown by Baxter (1994), twisted toroidal boundary 
conditions require little additional effort. 

2.2 Particular cases: Ising, Parafermionic 

For N = 2 the sum over I has only a single term. The angles <f>, ip drop out, 
Zj and Xj become Pauli matrices and we obtain the Ising quantum chain. If 
we put (j> = ip = then H describes the Fateev-Zamolodchikov parity invariant 
parafermionic Zjy-quantum chain. For N > 3 and <f> ^ or (p ^ parity is 
broken (therefore </> and ip are called chiral angles) and this gives rise to several 
interesting features of the model, e.g. to the appearance of incommensurate 
phases and oscillating correlation functions. 
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2.3 Yang-Baxter-integrable case 

If cos ip — A cos <p then H can be derived from a two-dimensional lattice model 
denned by the following transfer matrix which depends on rapidities p and q: 

L 

T P , q m, {i'}) = II w *Si w P Ah I's+i) (i) 

{/} and {/'} are the sets of iV-valued spin variables at alternating vertex rows of 
the diagonally drawn lattice. The rapidity lines p run horizontally on the dual 
lattice, the q vertically. The Boltzmann weights W and W are defined in terms 
of the functions x p , y p , ii v \ x q , y q , fj, q of the rapidities p and q by (Au-Yang et 
al. 1987) 

\ n n a n a 



M 9 / fJl y p -U J X q ' AX y g - uj3y p 

where n = 0, 1, N — 1. The requirement of Zjy-symmetry imposes several 
restrictions on the six functions, e.g. 

A(« + l) = < + <. 

where A is a parameter describing the inverse temperature. At fixed A only one 
of the three functions x q , y q , \i q is independent, the other two are determined 
by nonlinear relations. 

In the Ising case N = 2, the functions x q etc. can be simply expressed in 
terms of the meromorphic elliptic functions of modulus A: 



x q = — VA snq; y q — V\ enq/dnq, fj, g = yl — A 2 /dng 

and one obtains the parametrization of the Ising Boltzmann weights known e.g. 
from Baxter (1982). Integrability is guaranteed by the following Yang-Baxter- 
equation, which contains the p- and g-variables separately 

N-l 

W pA 1 ' - l)W r , q (l" ~ l)W r , p {l - I'") 



1=0 



The explicit expression for R pqr can be found e.g. in Albertini et al. (1989). The 
superintegrable case corresponds to fixing the p-dependent functions to x p = y p 
and Hp — 1 so that we are left with only one independent function, e.g. x q . 
H is obtained if at fixed p we expand the transfer matrix for small q — p. 
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2.4 The superintegrable case, Onsager's algebra 

The choice <j> = if = ir/2 gives the "superintegrable" case of the hamiltonian. 
Writing H = \N (A + XA t ) we have 

L N-l 7 i L N-l Y i vN-l 

3=1 1=1 j=l 1=1 

It has been shown by v. Gehlen and Rittenberg (1985) that A and A\ satisfy 
the Dolan and Grady (1982) conditions 

[Ao,[A ,[A ,Ai]}] = l6[A ,Ai}; [A u [A u [A u A Q ] } } = 16 [A u A Q ] . (3) 

Due to these conditions, A n and A\ generate Onsager's algebra A: 

[A u A m ] = IGt-m 

[Gi,A m ]=2A m+ i-2A m _ l ; [G l ,G m ]=0; l.mel. 

which has been an essential tool in the Onsager (1944) original solution of the 
Ising model. A implies an infinite set of constraints: 

[A 2 ,Ai] = [Ai,Ao} = [A ,A_i} = ... 

and the existence of the infinite set of commuting charges Q m : 

Qm = \ (A m + A-m + A(A m+ i + A_ m+ i)) with Qo = H. 

Finite dimensional representations of A are obtained if there is a recurrence 
relation of finite length among the A\ (Davies 1990, Roan 1991): 

n 

^ a k A k - n =0. 

k— — n 

If Zj (j = 1, . . . , tue) are the zeros of the polynomial f(z) = Y^k=-n °tkZ k+n -> 
then we can express the Am and G m in terms of a set of sZ(2, C)-generators 

n n 

A m = 2j2 (z>?E+ + zJ™Ej)- G m = J2 {*?-*J m ) H3 

where 

\ E ti E k] = S 3k H k; [Hj, E±] = ±25 jk E k t . 

So, A is isomorphic to a subalgebra of the loop algebra of a direct sum of sl(2, C) 
algebras. 

Each sector of the model is characterized by a set of points zj, which, of 
course, are not fixed by the algebra A alone. For the hamiltonian (2) the polyno- 
mials f(z), or rather related polynomials Vq(s), have been obtained by Albertini 
et al. (1989) and Baxter (1994) from functional equations for the transfer matrix 
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(1) of the two-dimensional model, as will be reviewed in the next section. 

As B. Davies (1990) has shown, from the property (3) of H it follows that all 

eigenvalues E of H depend on A in the form 

ros 

E = a + bX + Nj2±^l + 2\cos8 3 +\ 2 (4) 

3=1 

where a and b are integers and cos 0j ~ \{ z j + z J l ) ■ So the main task of the 
diagonalization of the superintegrable H reduces to the calculation of the cos 0j 
or equivalently, of the zeros of the polynomials f(z). For hermitian hamiltonians, 
the Zj must be on the unit circle. 



3 Solution of the Z 3 -superintegrable model 

3.1 Functional equation and ansatz of Albertini et al. 

Albertini et al. (1989) discovered the following functional equation for the Z 3 

( \ x / 6 

superintegrable case with the p- variables fixed to x p — y p = I y-^j J : Splitting 
off some poles from T Pj9 they get a meromorphic function T q which satisfies 

f q f Rq f mq = 3 L e~ iP { {t-l) L (uj 2 t-l) L f q 

+ (ujt - l) L (cu 2 t - l) L f R 2 q + {t- l) L (tut - 1) L f Riq } (5) 

where P is the momentum operator, t = x q y q /Xp and R is the mapping 

R(x q , y q , fig) = (y q , ujx q , fj,~ 1 ) . Solutions to this functional equation arc then 

obtained from the ansatz 

^= •■■ n (ttSt) n {d - a )(^+d± 2^(^-1)}. (6) 

i=i ^ ' j=i 

We skip here writing some factors relevant for the determination of the linear 
terms in (4) and chose a certain integer Pb = 0: for details see Albertini and 
McCoy (1990) or Baxter (1993,1994). The ansatz (6) contains mp excitations 
with Bethe-type rapidities vi, and in the last factor, mj functions Wj. It requires 
some calculation to see that if we put 




Wj = 

then the left hand side of (5) vanishes, and in order to get the right hand side 
vanishing too, tj must be zero of the polynomial 

j.—c ( m P i i j- m P 1 i J. 

1 I tJ 2 , . , -.xLTT 1 + v l t , , 2. , \L TT 1 + UV J^ 

t 3 



1=1 ■ ~ i ' 1=1 
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The integer c is chosen such that Vq becomes invariant against t — > u>t and Vq 
depends only on s = t 3 with Vq(0) finite. Despite its denominators, Vq{s) will 
be a polynomial if the vi satisfy the Bcthe-type equations: 

These Bethe-type equations resemble those of the spin-1 XXZ-chain. 
As mentioned already above, the hamiltonian is obtained by expanding the trans- 
fer matrix for small p — q around a fixed value of p. If Vq has degree mg, then 
from the tue zeros of (7) we get a set of 2 niE energy eigenvalues 

E t = o(l + A) + b - 6 (±tui ±w 2 ±...± w mE ) . (8) 

where e.g. for the lowest sector a = 3mE — 2L, b = 2Q. Comparing to the formula 
eq. (4) obtained from the Onsager algebra, we see that cos9j = (l+t 3 )/(l — t 3 ). 
Our hamiltonian is hermitian and so we should have | cos 9j \ < 1 and all zeros 
Sj = t 3 must be on the negative real s-axis. 

A difficult question is whether all solutions to these equations really give 
eigenvalues of our H, or whether there are spurious solutions. This problem has 
been investigated in detail by Dasmahapatra et al. (1993). 

3.2 Example: Sector m p — 

In general, in order to obtain energy eigenvalues of H for chain length L in the 
charge sector Q, one first has to choose a sector of given nip and to solve the 
Bethe-type equations obtaining a set of Bethe-rapidities vi . Then one has to use 
these vi to build up Vq and then to calculate its zeros Sj. 

Obviously, the simplest case is the sector m P = where there are no vi and 
one can directly start with building the Vq and looking for the zeros Sj. E.g. for 
Zj, with L — 20 and Q — 1 we calculate first 

Vi(s) = 20 s 13 + 8455 s 12 + 484500 s 11 + 8533660 s 10 + 61757600 s 9 
+ 210859245 s 8 + 363985680 s 7 + 326527350 s 6 + 151419816 s 5 
+ 34880770 s 4 + 3656360 s 3 + 146490 s 2 + 1520 s + 1, 

then solve numerically for its 13 zeros. Since the largest appearing coefficients 
grow exponentially with L (all coefficients are positive, and from (7) with mp = 
we get Vq(1) = 3 L_1 ), this procedure becomes impractical for L ^ 50. and e.g. 
energy eigenvalues for L = 1000 can hardly be obtained this way. Can we get an 
analytic expression for the zeros of the V{s) for large finite LI 

Baxter (1988) and Albertini et al. (1989) have obtained closed expressions 
for energy levels in the limit L — > oo. One method to get these is to rewrite the 
sum over the terms determined by the zeros (chosing all ±-signs to be + and so 
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obtaining the lowest eigenvalue E) as an integral over a contour C enclosing the 
negative real t-axis: 



E \/ 1 + 2Acos ^ + A2 = X dt V (1 - A)2 + — s ln7W 

Then the contour is opened and gets caught by three 120° symmetrically oriented 
cuts, which arise from the square root of the integrand. The rightmost of these 
cuts is on the positive real axis at 1 < t < |(1 + A)/(l — A)! 1 / 3 . All three 
cuts give the same contribution, and the rightmost one can be evaluated easily, 
since for L — ► oo only the term proportional (t 2 + t + 1) L of V, eq.(7), is 
non-oscillating and needs to be kept. Baxter (1988) finds that the ground-state 
energy for general Z N is given by (see also Albertini et al. 1989): 



,. E N /■^ 1 " JV > n sin™ (x + § )-sin"x x9 

lim — = / dx*l-2X ^ ^ s^ + A 2 . 9 

4 Determination of the zeros of the 'Pq(s) 

We attempt to find the zeros of the Tq(s) for finite L on the negative real s-axis 
analytically. From the explicit expressions and numerical studies we observe: 

• For L = Q mod 3 the zeros come in pairs which are reflection symmetric with 
respect to s — — 1: with Sj also 1/sj is a zero. This remains approximately 
(within a shift) true for the other charge sectors. 

• The distances between adjacent zeros alternate slightly in size so that a 
smooth curve Sj(k) is obtained only if we connect every second zero. 

4.1 Case Z 3 , m P — 

A useful variable transformation to be applied to Vq (s) turns out to be 
t = (u— \/3) /(u + VS) or, writing u — tan [3: 

t _ sinQ9-f) 



sin (/?+!)• 

This maps -oo < s < to e.g. -f < (3 < § or ^ < [3 < *f (for 
definiteness let us choose the first interval), giving 

9 3 , 9 —3lu cos/?e~ I/3 

t 2 +t + l = Ty and t A +ut + uj= „-. 

(2sin(/3+f)) 2 2(sin(/3+f)) 2 

so that 

V Q ~(t 2 +t+ l) L + UjQ(t 2 + UJ 2 t + L0) L + LO-Q(t 2 +LOt + L0 2 ) L = 0, 

multiplying by the non- vanishing factors (sin(/3+ ^)) becomes 
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L= 1 2 



Q=0, 1 




tt/3 



Fig. 1. Left-hand (dash-dotted line) and right hand (dotted and dashed) parts of (10) 
for L = 12 sites. Charge sector Q = 0: dotted lines, sector Q = 1: dashed lines. The 
intersections with the dash-dotted line with < f3 < 7r/3 give the me solutions 

Pk- ttie = 8 for Q = and the = 7 for Q = 1. 



-(2cos/3)- L = 2cos(l{/?+|} + ^). (10) 

Now - 1 < /3 < | means 2 cos /3 > 1 so that for L » 1 the left hand side 
of (10) will be small, except near the corners of the interval in /3, as is shown 
in Fig. 1. In order to get a solution, the fast oscillating right hand term has to 
coincide with the non oscillating left hand term. We see that all but the solutions 
near (3k ~ — f and (3^ ~ § are practically equally spaced in the variable (3. 
Neglecting the left-hand side altogether, we get the approximate solutions 

h^K-*-- K= 6k+ ^~\ ; k = l,2,...,m E . (11) 



Finite-size energy levels of the superintegrable chiral Potts model 



9 



or, expressed in terms of the cos 6^ which we need for the energy eigenvalues: 

1 + 4 sin 3 (^+f)-sin 3 ^ 
cos# fe = -;—£ = -_ (12) 



sin 3 (if + §) + sm 6 K 



It is easy to see that using (11, 12) and replacing the summation in (8) by an 
integral leads to Baxter's formula (9) for the thermodynamic limit. Through a 
Q-dependent shift (11) describes most of the Q-dependence of the finite-!/ levels. 
Further finite-size effects arc due to the deviations of the largest and smallest 
(3k from the approximation (11). 

In order get an idea of the precision of (11), consider e.g. L = 30, Q = 0. With 
d k = {(3 a k ppr / P e k xact ) ~ 1 where PT Pr is the approximation (11), we find 

di = 0.0012; d 2 = -1.89 • 10" 5 ; d 3 = 6.48 • 10~ 7 ; d 4 = -4.39 • 10~ 8 ; . . . 

For L = 30000 and Q = we get d x = 0.99 • 10~ 6 , but compared to L = 30 also 
the spacing between neighbouring (3k is ~ 2 • 1000/3 times smaller. 

Taking the arcsin of (10) we get a pair of useful (still exact) equations for 
the (3 k : 



=F arcsin ■ 



2ttQ 
3V 3 



^^2 



2 (2 cos (3 k ) L 

The two sets of counting integers I k vary within the limits 



0. 



(13) 



with 



3Tl 



It = 



L+Q-lTl 



(14) 



The arcsin-equation has the great advantage that from it each zero can be cal- 
culated separately by fixing Numerically on a small PC we can easily get 
zeros with 40 digits precision for up to L = 300000 sites. Recall that comput- 
ing directly the zeros of the "Pg(i) we were limited to L <, 50. In addition, 
the arcsin-equation lends itself to excellent analytical approximations for the 
"corner" -solutions: Let us call £ the term which is important only close to the 
corners: 

^ arCSin 2(2cos/^ (15) 

We find numerically that £^ -C 1 even for the extreme cases = I^b- We use 
this feature to approximate . Let us consider the lower end of the interval in (3, 
i.e. I k I a - Then we get a behaviour smooth in L when keeping Q fixed (At the 
upper end I k ^ h the behaviour is smooth only if we keep L + Q fixed mod 3). 
Using arcsinx w x and writing (3 = — n/3 + e we get (omitting the superscripts 



(2 cos/3)~ L = \ (cos e + V3 sine) w \ (l + VTh 



-L 



(16) 
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On the other hand, solving (13) for the (3k which is not in £fc, we have: 

£fe = £ Wlth fc = I fc ~ ~3~ 2 / ^ ' 

Inserting (17) into (16) and assuming L> 1 we get 

& « ^ (l + 6) ) ^ 5 exp {-V3(A ± 6)} (18) 

This nonlinear equation for valid for small e, can be solved by iteration. The 
zeroth and first iterations give 

4 ±(0) - \ and = \ exp {-V~3 (if ± \ e"^) }, 

respectively, or 

So we find corrections to (11) which fast become smaller with increasing 7fe, and 
are useful mainly for two smallest values (3k ■ Closer to the middle of the interval, 
anyway the approximation = gets very accurate. 

Coming back to the example L = 30000, Q = given above, adding the ith 
iteration corrections /L improves the approximation (11) to 

df ] = -5.48 • 10~ 8 , = 3.285 • 10" 9 , <i^ 2) = 1.53 • lO" 10 
for the zeroth, first and second iteration, respectively. Regarding the (3k away 
from the corners, already from (11) we get e.g. d 7 = 1.4 • 10~ 25 , which improves 
slightly to (iy -* = 6.3 • 10~ 26 by adding the zeroth iteration. 

Let us study how much the results for the are affected by the approxima- 
tions made in (16) and (18), and how many iterations are needed. In Table 1 
we give numerical values from our iteration formula in order to compare them 
to the following numerical results from the exact formula (13). A fit of the In- 
dependence of the numerical solution, obtained from L = 25000, . . . , 125000 
sites is, expressed in terms C ± / 7r: 



Q = 0+ 
= 0- 
= 1+ 
= 1- 

= 2+ 
= 2- 



+0.009927623657298 + 0.04835895/L + 0.0315/L 2 + . . . 
-0.000045412321296 - 0.002017284/L - 0.0338/L 2 
+0.00000740160385053 + 0.0004910486/L + 0.0130/L 2 + . 
-0.0017242062241336 - 0.023760499/L - 0.0922/L 2 
+0.0002780560828110 + 0.007462876/L + 0.0685/L 2 + ... 
-0.00000120676616449 - 0.000111825/L - 0.00430/L 



2 



We see that for large L our equation (18) determines the correction (15) to (11) 
with excellent precision if its twice iterated solution is used. The coefficients in 
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these expansions indicate that the asymptotic large- L- behaviour sets in only for 
about L k, 1000. 

If we move downwards from the upper corner = tt/3, we get exactly the 
same numbers and formulae, just signs are changed and charge sectors are per- 
muted. 



Table 1. Values of £k for the lowest value of k , as obtained from (18) at different 
iteration levels. We have multiplied the values by 10 6 /n in order to avoid too small 
numbers. 



Iteration 


Q=0+ Q = 0- 


Q=l+ Q=l- 


Q=2+ Q = 2- 


Ost 
1st 
2nd 


10477.0 -45.401101 
9896.39 -45.412318 
9927.70 -45.412321 


7.401901956 -1708.097 
7.401603838 -1724.047 
7.401603850 -1724.196 


278.47707 -1.20675824029 
278.05541 -1.20676616444 
278.05605 -1.20676616449 



4.2 The leading finite-size corrections 

The formulae of the last section are sufficient to calculate finite-size eigenvalues. 
But often it is convenient to transform the sums in (8) into integrals in order 
to have a compact expression also for finite L eigenvalues or to get the lead- 
ing finite-size corrections to the thermodynamic limit results. This can be done 
e.g. applying Euler-Maclaurin-tcchniqucs as is common in calculations for the 
XXZ-quantum chain. Following Woynarovich and Eckle (1987) and Hamer et al. 
(1987), but taking care of the alternating spacing of the zeros, we define two 
functions Z^ 

zm . ^{ Tarcsi »_L_ + £ (, + f) + 2f» 

where — 3m^/L so that hin^^oo^ = 1 . From (14) we read off m E = 
if — if + 1. In Zf the roots 0k are exactly equidistant. We define the root 
density 

and use cr+ (0) + <r^ (/3) = CTqo, mf, + m E = the- Then the finite-size 
correction to a quantity J d0 w{0) (since we are mainly interested in the 
ground-state energy per site, this will be the w of (8)) is 

( m E m E \ 7r/3 

E H#) + £H9fc) - / d0a oo w(0) 
k=l kTl / J-*/* 
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w(J3) 



/■tt/3 / rOt fir/3\ 

= dPiat+aZ-a^wiP) / +/ + ^ CT 

+ ^ ?Wft±) + ^ )H _L ? (^_^) + ... (19) 

Both /?± = /3f = /3±(/±) and ^ = P ± ( I b) depend on the charge g. 

5 The zeros for higher Z^-superintegrable cases 

Analogous trigonometic formulae result for the higher Zjy-cases. With increasing 
N, the equations get more involved, but have a similar structure. 



5.1 The case Z 4 : 

Since we transform the variable t and not s = t 4 , in order to map the negative 
real s-axis to a range of real (3, we have to include a phase and use 

t = Sin j^'|; (20) 
sin(/3+f) V ' 

so that — oo < s < is mapped to e.g. — j < (3 < j. After some algebra, we get 
that the condition for the zeros of the mp — 0- functions Vq(s) 

Pq(s) = t~° {(t 3 + t 2 +t+ l) L + (-l) Q (t 3 -t 2 + t-l) L 

+r Q (t 3 + it 2 -t- if + i Q (t 3 -it 2 -t + i) L ) = o (21) 

translates to 

(5^r-H- + i) + ?} + -{<-^) +! ?} 

= (22) 

where p± = arctan { (-\/2 ± l) tan (3} . 



For L >• 1 and cos (2/3) not too close to zero, the second term of (22) dominates 
and the solutions come where 

cos{L(p + + f) + ^2} « 0. 
After some rearrangement this leads to 

tan (3 k ~ - tan - tan - cot ( if + - with if = -f- — , (23) 

4 8 V 8 / AL 

k = 1,2,..., [(3L — Q)/4] and e.g. — j < (3k < f. In contrast to the Z 3 -case, 
where in the approximation (11) the zeros were equally spaced in (3, here this is 
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no longer the case: now the f3k are more dense in the center region. In the center 
(3 w the oscillations of the p+-term in (22) are (V2 + l)/(\/2— 1) ~ 5.83 times 
faster than those of the p_-term. (20) and (23) are equivalent to 

1+4 sin 4 (if + f ) - sin 4 (if ) 

" ~ 1-3 " sin 4 (if + f)+sin 4 (if)- 

As in the Z 3 -case, in the thermodynamic limit, this approximation leads back 
to Baxter's ground-state energy formula. 



5.2 The case Z 5 : 

Transforming t = sin (/? — ^)/ sin (/? + J), the condition 

Tq ~ J L + u- Q jf + w Q j£ + u;- 2 Qj 2 L + lu 2 ®J 3 l = 
with J k =t 4 + cj k t 3 + uj 2k t 2 + Lu 3k t + Lu 4k takes the form 

1 /C0S/?\~ L f / 7T\ 27T<3 

+ cos { L 



2 Vcosf 







+ ^ 2 cos{l( X+ |) + ^}.0. (24) 

7T 7T n , „ 4sin 2 /3 + V5 - 2 

where tan y = cot — cot — tan p and B — . 

5 10 2 cos | 

The solutions come for — oo < s < corresponding to e.g. — j < (3 < j. Since 
for L > 1 and not too close to f we have 1 < (cos/3/ cos |)~ L < B~ L / 2 , 
the solutions come close to where the cosine containing \ vanishes, i.e. at 

ir ir n ^ (lOfc + 2Q - 5)tt 
tan/3 fc w - tan - tan — cot \K + —) ; K = — 

and integers k = 1, 2, . . . such that e.g. |/3fc| < -|. Equivalcntly, we can write 



mSk . 1+4 „ .^^1-*;' (25) 

sin 5 (if +f) + sin 5 if 



The i^-formula (10) has the same form as the first terms of (24) since 2cos/3 = 
cos (3/ cos j. The expressions for general Zn will be given elsewhere. 
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5.3 The case Z 6 : 

We conclude with the case Z 6 : As in the case of Z4, because N = 6 is even, we 
include a phase and transform —00 < s < to — f < /3 < f using 

k ,w« sin (/? — £) 

sin 09+ f) 

Now Vq{s) = (we don't give the explicit Z§ formula for "Pq, since this is 
easily generalized from the Z 5 -case) after some algebra leads to 

+ s) + ¥} + -^M 1 (« + 5) + f } 

+ coJl(W+SV^) = (26) 



V3 tan 2 J - tan 2 /3 

with a = — ~% g— and 

2 tan 2 f+tan 2 /3 

po = arctan^V3 tan/?); p± = arctan | (2 ± V3~) V3 tan/3 j. (27) 

For L 3> 1 and a not too close to zero, the last term containing Lp + 
dominates and the solutions (3k approach 

- 7T 7T / 7T \ (6fc + Q — 3)7T 

tan/3fc w -tan -tan — cot + — J; if = — 

and k = 1,2, .. . ,die with raj = [(5i — Q)/6}. We can rewrite this like (25). 



6 Conclusions 

In the analytic solution of the superintegrable chiral Zjv-Potts quantum chain 
the energy eigenvalues are determined by the zeros of special polynomials. These 
polynomials also play a central role in the representation theory of the Onsager 
algebra. In the thermodynamic limit, sums over functions of the zeros can be 
evaluated by contour integrations. For finite systems we can try to calculate the 
zeros numerically. However, applying standard computer routines directly to the 
explicit form of the polynomials will run into trouble beyond some chain length 
L, because the largest coefficients of the polynomials are growing exponentially 
with increasing L. So we were looking for an alternative and possibly analytic 
approach. 

By a suitable change of variable we transform the equations for the zeros 
into trigonometric equations, for which approximate analytic solutions can be 
obtained. We give explicit formulae for the simplest polynomials. 



Finite-size energy levels of the superintegrable chiral Potts model 
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In the case of i? 3 the solutions corresponding to the zeros of the polynomials arc 
almost equidistant in the variable (3. However, for Z 4 and higher, the density of 
zeros in f3 varies appreciably over the finite interval. Our formula describes the 
solutions, except those near at the boundaries, with exponential precision in L. 
The leading finite-size effects are due to the solutions at the extreme ends of the 
interval in f3. For the case Z 3 we derive a special formula which describes also 
the extreme solutions, thus allowing to calculate finite-size energy eigenvalues 
for the low-lying sectors. 
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